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<N ! Abstract 

> ■ 

We describe the systematical construction of the first order formalism for multi-centered 
\^ ■ black holes with flat three dimensional base-space, within the so-called T 3 model of N = 2, 

If) . D = 4 ungauged Maxwell-Einstein supergravity. 

I/-) I The three possible flow classes (BPS, composite non-BPS and almost BPS) are analyzed 

in detail, and various solutions, such as single-centered (static or under-rotating) and all 
known multi-centered black holes, are recovered in this unified framework. We also consider 
the possibility of obtaining new solutions. 

• 1— 1 ' The almost BPS class is proved to split into two general sub-classes, corresponding to a 

positive or negative value of the duality-invariant polynomial for the total charge; the well 
known almost BPS system is shown to be a particular solution of the second sub-class. 



1 Introduction 



Multi-centered non supersymmetric extremal black holes have been the subject of a number of 
recent studies [3]- [8], one possible reason being that supersymmetry actually plays a less crucial 
role than extremality in the investigation of these solutions. 

The supersymmetric class for multi-centered black holes was investigated by Denef [9], in 
order to match the true BPS spectrum and the spectrum of spherically symmetric black holes in 
supergravity. The resulting multi-centered BPS first order equations, whose solutions and prop- 
erties were then studied in [10], acquire the well-known form after integrating out the auxiliary 
phase field in the single-centered limit. 

More than ten years passed before a systematical construction of solutions for the under- 
rotating non-BPS branch was worked out by exploiting group-theoretical techniques [3, 4]. 

The present paper is devoted to the determination of the most general flow equations under- 
lying the first order formalism [11]- [13] for non-BPS multi-centered and/or rotating black holes 
with flat three dimensional base-space. For simplicity's sake, we will consider the so-called T 3 
model of N = 2, D = 4 ungauged Maxwell-Einstein supergravity. But we expect all the results 
to be generalizable to the STU model, a task left to further forthcoming investigations. 

The first attempts to write down flow equations for non-supersymmetric multi-centered con- 
figurations were made in [8]. The method was based on the squaring action procedure, but 
unfortunately an explicit form of flow equations in terms of charges, fields and other parameters 
was missing. In the present investigation, we try to fill this gap, and write down first order 
equations for non-supersymmetric multi-centered configurations in an explicit form. 

Besides its inner simplicity and elegance, the importance first order formalism relies in the 
possibility to switch from second order differential equations of motion to the first order ones, 
without doubling their number. Indeed, even if some auxiliary fields are added in the multi- 
centered case, they turn out to satisfy algebraic equations. This is due to the fact that the scalar 
charges are not independent, and the formalism automatically discards the blowing up solutions. 
Integrating first order differential equations is certainly easier, and thus the possibilities to find 
the corresponding attractor flows are considerably enhanced. In this paper, we will retrieve a 
number of known solutions in a unifying setting, and we will predict the existence of other new 
ones, at least at the level of the most general system of equations governing the flow. 

The paper is organized as follows. 

The setup and notations are defined in Sec. 2. 

Then, Sec. 3 presents the general method of construction of first order equations, considering 
the well known BPS class. 

In Sec. 4, this method is then applied to the construction of the first order formalism for the 
composite non-BPS class, and well known examples, such as single-centered static and under- 
rotating black holes, and multi-centered electric and magnetic solutions, are recovered. 

Sec. 5 exploits the same procedure for the broadest class of almost BPS mult i- centered 
solutions. The existence of two general sub-classes of first order flow equations is proved, and 
the well known almost BPS system [5] is retrieved as a particular solution of the second class. 

Finally, Sec. 6 contains a summary and an outlook of results. 
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2 Setup and Notations 



The bosonic part of the iV = 2, D = 4 Maxwell-Einstein supergravity Lagrangian density action 
for the so-called T 3 model reads 

= -\ R + <f> G djdj+ \ tMAsF^F^ + \ u^F A * F s ^, 



where F A = dA A is the electromagnetic field, G = — uzm is the metric on special Kahler 
symmetric scalar manifold SL(2, M)/U(l) , and the scalar- vector couplings are given by 



i 



(t-t){t 2 + Ut + t 1 ) -3(t 2 -?)\ 1 { ~(t + i) 3 3(t + t) 
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~ 4 V — 3(t 2 — t 2 ) 6(t-t) /' ^ 4 V 3(t + 1) 2 -12(t + t) 

The Ansatz of under- rotating multi-centered black hole (with flat D = 3 spatial slices), 
independently on the branches, reads 

ds 2 = g^dx^dx" = e 2U (dt + uM) 2 - e- 2U 5 ij dx i dx j , (1) 

where % denotes the Cartesian indices, and /i, v are space-time indices. The functions U, u>i and 
the complex modulus field t = x — iy, where x is axion and y is dilaton, are solutions of the 
following equations [1]: 

dift = 0, & d t [e~ 2U M a p [djiP - fjbP)] = 0, 
Gd {i td j) t + diUdjU + ± e 4U fifj = - 1 - e-WdiVMapdjlP, 
8" dftU + l - ' "'d' 1 /;/, = - 1 - e-^d^M^djbP, (2) 



8 ij di (Gdd) - ^ S^ditdd = - e- 2U 5 ij dib a ^^ 9,, 
v ' dt 3 2 dt 3 



where a, (3 are symplectic indices. 

In the present framework, /, corresponds to the rotation, and it is defined as follows: 

/. = S^dkUL (3) 

On the other hand, b a denotes the electromagnetic potential. As usual, we define the dual field 
strength Ga = dC/dF A and the whole field strengths' symplectic vector F a = [F A , Ga), so that 
the electric component of the field strength F a reads 

K = ~dib a . (4) 

As for M a p, it is a positive definite symplectic symmetric matrix, which defines the single-centered 
black hole potential [2] and reads 



j"AS + ^AA'At A ' E '^S'E ^AA'Ai A ' E 
,, AA ',/ ,, AE 
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where ^ae^ A ' = ^a'- 

A crucial point is the integration of the electromagnetic field equations. In the following 
treatment, we will use the following Ansatz: 

e- 2U M a pd^ = n aP Ql (5) 
where Q? = (d^ + 0) and Q a p is the bilinear symplectic invariant structure: 

-6f 



n 



a/3 



'A 



^ o 



and H 13 are harmonic functions containing only appropriate electric and magnetic black hole 
charges (namely, fluxes of the 2- form field strengths). 

It is worth remarking here that one may generally write Cf { = (h? + filrj , where iff is 

the dual of the three dimensional field strength (not necessarily harmonic), and S^djH? = 0. 
However, we will see below that (5) is effective for all classes of multi-centered black holes and 
models under consideration. 

Substituting d^P given by (5) back into equations (2), one obtains 

F difj = 0, e- 2U M a(S d,^ = n aP Ql 

1 
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G fytfyt + diUdjU + - A e w fj, = e 2C %, 



S ij <),<) s r + l - e AU 5^fJ 3 = <-' <S ; * \) r (6) 



v ' dt 3 dt 



where 

Vij = Z(iZj) + Z t (i Ztj), (7) 

is a generalized black hole potential, and (v a denotes the covariantly holomorphic section of 
special Kahler geometry throughout) 

Zi = Qf Otf v\ Z ti = E l D t Zi = Q? Q aP vl, (8) 

are Cartesian-indexed generalization of the N = 2 cental charge function and of its first Kahler- 
covariant derivative 

Z = Q a n aP v p , Z t = E l D t Z = Q a tt aP v? , (9) 

in which the charge vector Q a = {p A , q\} is replaced by Qf . E l = £| is the (inverse) Vielbein 

of the scalar manifold (\E t \ 2 = G); and we henceforth refrain from hatting the subscript "t", 
denoting the flattened scalar coordinate. 

Is worth mentioning here that an object like the generalized black hole potential Vij appeared 
in [14], but our definition has a more transparent geometrical interpretation. As for single center 
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static and spherically symmetric black holes, only the V TT (where r = 1/r is the inverse radial 
coordinate) component survives in this case, and it equals the usual black hole potential Vbh [2]. 

We are now going to construct the first order flows for these equations in various branches. 
After [3], it is known that for under-rotating multi-centered black holes three possible flow classes 
exist, namely: BPS, non-BPS composite and almost BPS. 

For all classes of multi-centered flows based on the known solutions [3], [4] and [9], one can 
consider the following Ansatz for the first order equations: 

diU = mfe, z u , i ;,«„) e u , (10) 

dit = tffiifZi, Z u , fi,a a }e u , (11) 

diOl a = T ia fe, Z ti , fi,0£a) ^ ■ (12) 

Here Wi denotes a Cartesian-indexed generalization of the first order (fake) superpotential W 
for single center black hole [11]. flj is the corresponding generalization of the scalar charges 
(2E t dfW), and en's are some auxiliary fields. Qa va sans dire, these functions will generally be 
different for the various classes. 

For pedagogical reasons, we will start by considering the BPS class. 



3 BPS Flow 

This is a well known case [9], and we will consider it in order to establish the procedure which 
we will exploit for the non-supersymmetric cases. 

For single-centered BPS black holes in spherical coordinates, the first order superpotential 
reads 

W = \fz~2. (13) 

In order to re-express it in Cartesian coordinates, a mere replacement of charges by derivative of 
harmonic functions (namely, Z — > Zj) is not convenient, because Cartesian indices would occur 
under square root, and covariance would be lost. Actually, manifest covariance can be preserved 
by introducing the phase a and expressing W (13) as [9] 



W = Re [e~ iao Z] , 




Then, one can easily switch to Cartesian coordinates by introducing the corresponding index: 

Wi = Re [e~ iao Zi] . (15) 

Of course, the definition of «o through Z and Z does not hold true anymore here, and the 
corresponding equation should be established. Rotation can be introduced into by prolonging 
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central charge as Z{ — > Zi (recall (8)). By exploiting such a procedure also for the scalar charges, 
one ends up with the following first order system: 



diU = Re e~ iao Zi e u , (16) 
dit = E'e^Zue 17 . (17) 
By plugging these equations back into (6), one then gets equations for a and fi, namely: 



dia 
fi 



e u Im 
2Irn 



-l ao 



Vsz t 



-ia 



(18) 
(19) 



We notice that Eq. (19) contains f\ in both sides; thus, by recalling definition (8), it can be 
rewritten as follows: 



1 - 21m [e~ ia ° b a Vt aP v p ] e~ u )fi = 21m [e~ ia ° Zj e 



(20) 



It is here worth noting that, due to the presence of the rotation, all equations contain the 
electromagnetic potential. One can easily show that the following Ansatz solves the equation for 
the electromagnetic potential in a mo del- independent way: 



b p = 2Re [e~ iao v?] e u . 



(21) 



We point out that this Ansatz can be constructed by replacing all Zj's with the symplectic 
sections v 13 in the superpotential Wi, and then by multiplying by 2e u . In the following sections, 
we will see this procedure to (partially) work also within non-supersymmetric classes of flows. 
Plugging (21) back into all other equations, one finally gets the Denef system [9]: 



diU = 


Re [e 


~ iao Zi] e 




= 


E l e iao Z u e h \ 




dia = 


— Im 


dt 9 * 


+ 


fi = 


-21m [e~ iao Z 


i] e 



2 f '' 
-u 



„2U 



(22) 
(23) 

(24) 
(25) 



where K = — log (—i(t — t) 3 ) is the Kahler potential. 

The rotation should be divergence-free (namely, 5 l i difj = 0), thus so from (25) one obtains 



fi = diH a VL a s . 



(26) 



Let us now consider the static single-centered limit of the system (22)-(25). By setting fi = 
from the start, the standard BPS solution is obtained without any integration of the equation of 
electromagnetic fields; however, the condition 

Q a n a ph,P = 0, (27) 
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where h 13 are the asymptotical values of the harmonic functions H@, should be introduced ad 
hoc. 

On the other hand, after integrating electromagnetic fields out and finding (26), by setting 
fi = one obtains the condition (27) automatically. 



4 Composite non-BPS Flow 

4.1 General Analysis 

In this section we are going to investigate the so-called composite non-BPS branch [3]. In order 
to determine the fake super potential and the corresponding scalar charges, in this case one can 
start from the single-centered non-BPS black hole fake superpotential [13], and introduce one or 
more phases. Due to the rather cumbersome explicit expression of the non-BPS W, it is not easy 
to guess how to consistently introduce phases, so we choose to proceed in another way. Namely, 
we will use the procedure described for the N = 2, D = 4 STU model in [3]. As treated in that 
paper, under-rotating solutions correspond to various nilpotent orbits, and thus one can try to 
switch from the BPS orbit to the non-BPS ones by performing the following steps (for the T 3 
model): 

1. Consider the scalar charges and the quantity whose real part is the BPS superpotential 
W = Re(tw), and remove the phase cvq therein: 

w = Z, n = Z { . (28) 

2. Switch to the duality frame (on spatial infinity) in which the axion x vanishes and the 
dilaton y is set to 1, thus establishing a one to one correspondence between the the charges 
and the central charge and its Kahler-covariant derivative. In such a way, all quantities 
can be rewritten solely in terms of the charges themselves. 

3. For the composite non-BPS orbit, the signs of go and of the combination p° + qi/3 should 
be flipped, and such quantities should be written in terms of the central charge and of its 
derivative, thus yielding: 



w 



1 -(Z-V3Z-^ U = ~(V3Z + Z t y (29) 



4. Then, the phases should be introduced using the following rules under [[/(l)] 2 transforma- 
tions [3]: 

Z -> Ze^ {ao+3ai \ Z t ->■ Z t e~^ ao+ai \ w ^ w e~^ ao+3ai \ U^Ue^ a °- ai \ (30) 
thus achieving the expressions for the fake superpotential and corresponding scalar charges: 



W = -Re 
2 



e~ 3iai Z - V3e~ iai Z- t 



(31) 



n = ~(V3e iai Z + e~ iai Z t y (32) 
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Clearly, the expressions (31) and (32) are written only at spatial infinity. However, we 
choose to use them all along the flow, by just replacing the charges Q a with the quantity Qf, 
exhibiting the same symplectic covariance.This replacement can be understood by observing that 

naturally occurs in the second order equations (6), and at infinity it equals Q a when spherical 
symmetry is restored. After integrating the phase a\ out, we will show below that the known 
expression of non-BPS W for the T 3 model [13] is retrieved. 

By using the expressions of Wi and IT,, which as in the BPS class contain only one phase (in 
this case denoted as a\), the following first order equations can be written down: 



-ICtl 



Z t 



jj 



diU = -Re 
Bit = —E 1 (Vle^Z, 
Thence, by plugging them back into the second order equations (6), one obtains 

3 iQl EL 



diCti = -=e u sin oi\ Re 

v3 



fi 



Im 



-3ia 



1 Zi - V3e~ iai Z, 



-v 



(33) 
(34) 

(35) 
(36) 



It is here worth stressing that, as in the BPS class, fi is the imaginary part of the quantity whose 
real part is Wi. 

Moreover, /j enters both sides of (36); thus, we would like to integrate out the electromagnetic 
potential completely, as for the BPS class. Following the procedure outlined in Sec. 3, we replace 
the central charge and its flat derivative respectively by v a and vf, obtaining the following Ansatz: 



ti 3 = Re 



-iai 



(37) 



Unfortunately, this Ansatz does not work here well. Indeed, by inserting it into (36), it turns 
out that fi drops out, and one ends up with the expression: 



Im 



-3ia 



Z t - V3e~^ Z u 



0. 



(3* 



However, we can add to b° (37) a new term proportional to a boundary function b. As we will 
see, b represents a real rotation, whereas in the BPS class only intrinsic rotation fi is present. 
Consequently, the new Ansatz reads 



b^ = Re 



e -3i Qi v 



e u -be 3U Im 



-iai Trfi 



(39) 



The new term added, proportional to b, is nothing but fi (36) (namely, the imaginary part of 
the quantity whose real part is Wi) with the central charge and its flat derivative respectively 
replaced by and . Using this Ansatz, one finds that (38) holds true, whereas the equation 
for electromagnetic fields yields the following expression for ff. 



fi = dib + 26-^Im [e~ 3iai Zi] + V3be u Re V3e 3iai Z { 



D iai 



z 



t i 



(40) 
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Finally, by inserting (38) and (39) back into the other equations, one obtains full-fledged first 
order formalism for the composite non-BPS class of under-rotating multi-centered flows: 



d t U = -Re 
2 



e -3fai z . _ ^-im z u 



e U + \bf,e AU (41) 



dit = - l -e u - 2 ^E l (V3 e u (i - be 2U ) + v^e 3 "* 1 Z t + e iQl Z ti ) . (42) 
4.2 Some Examples 

In order to find explicit solutions, the system (41)-(42) should be supplemented by the condition 
(38), which is the integrated version of equation (35). 

1. Let us consider the single-centered non- rotating non-sup ersymmetic black hole in the T 3 
model. By setting /, = before the second integration of the electromagnetic field equations 
and bearing in mind that, by virtue of spherical symmetry, all three Cartesian equations 
in (38) can be replaced by the only radial coordinate, one obtains the following solution 
for ol\\ 

e2iai 3\ + t v ^9~X 2 + 2AXZ t Z- t + 16Z t Z- t (3ZZ-Z t Zt) | Z t 

4V3ZZt V3Z' 

t - 3 \ V 3 

ZZ - ^\ - h [ZZ + ZtZf) + -^j=i {ZZf - ZZf) 
ZZ - - h {ZZ + Z t Zt) - (ZZf - ZZf) ) , (44) 



1/3 



where I4 is quartic invariant. Plugging (43) back into multi-centered fake superpotential 
(31), the known expression for the non-BPS fake superpotential of the T 3 model [13] is 
retrieved: 

W 2 = ](ZZ + Z t Z t )+h. (45) 

Of course, one can also consider the single-centered limit after the second integration of 
electromagnetic potentials. In this case, one has to impose the condition fa — 0, where 
fi is defined by expression (40). By setting b to a constant, and integrating ot\ out, one 
achieves the following condition: 

(x 2 + y 2 ) (-1 + v 2 ) -a 2 _ + v 2 a\ + 2be 2u yu(a^ + a + ) - 2x + u 2 a+) = 0, (46) 

where x = Re (t) is the axion field, y = — Im (t) is the dilaton field, and cr + , a_ and 
v are defined in formula (3.19) of [13]. It is easy to check that solutions (4.19) of [13] 
indeed satisfy this condition. As for the BPS case, one obtains the condition (46) in 
order to have no rotation, whereas in the previous approach (namely, before integrating 
electromagnetic potential out) such condition had to be imposed by hand. Thus, one can 
conclude that the composite non-BPS class exhibits a smooth single-centered non-rotating 
non-BPS (extremal) limit. 
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2. The example concerns the generalization of Rasheed-Larsen black hole [15], namely a non- 
BPS single-centered rotating black hole with all charges switched on. As we are considering 
a single-centered black hole, all harmonic functions sourced by the charges depend only on 
the radial coordinate r. Thus, condition (38) can be solved as in the case of a non-rotating 
single-center black hole. By computing a% and inserting it back into the expressions for 
the derivatives of warp factor and for scalar fields, one obtains the very same equations for 
the non- rotating single-centered black hole, and the very same solutions (see e.g. (4.19) of 
[13]), the unique difference being that, b now is not constant anymore, but instead it is a 
harmonic function which represents rotation: 

b = J» + J^, (47) 

where bo, J are constants and J is the angular momentum. By setting all charges to 
zero but the DO and D6 ones (which, in the notation of [13] amounts to setting v — > 0, 
(j_ — )• and v<j + — > (qo) 1 ^ 3 / (p ) 1 ^ 3 ), and by imposing the additional condition hi = h 
(corresponding to the vanishing of the axion), the Rasheed-Larsen solution is consistently 
reproduced. 

3. As third example, we consider the ( T 3 -degeneration of the) composite multi-centered black 
holes investigated by Bossard and Ruef in [3], [4]. In the multi-centered case, the condition 
(38) can explicitly be solved only in some particular charge configurations. Namely, when 
all D6 charges or all DO charges vanish on all centers (respectively called magnetic and 
electric configurations). For the magnetic configuration, condition (38) yields «i = 0. By 
solving all other equations, one gets the following solutions for warp factor, axion and 
dilaton, respectively: 

U = -ilog[4^i(-i/ p i) 3 -6 2 ], (48) 
X = (49) 

-2U 

y = w (50) 

Here H p i < denotes the harmonic function which corresponds to the magnetic charges 
p 1 , while the harmonic function B satisfies 

di B = (51) 
H p i 

Furthermore, the non-harmonic function A is constrained to satisfy the equation 

diA = diH qo + 3 A + B d t H qi - 3 B 2 . 

It is here worth mentioning that one could have found this solution more easily by con- 
sidering equations for the combinations e~ u ' j ' ^fy and x — ybe 2U . As for fi, it reads 

f^dib+iHpifdiB, (52) 
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and b is such that implies fi to be divergence-free. For the electric configuration, condition 
(38) yields ot\ = arccot (x/y) and the following solutions can be obtained: 



U 



x 



-4 log 



27 



A(H, 



91 



f-b 2 



b-2/9(H qi ) 2 B 



2(-bB + (((H qi )^)/9 + (H Ql A)/3) 



2(-bB + (((H qi )^)/9+(H qi A)/3y 
where functions B (harmonic) and A (non-harmonic) satisfy: 

2 diH qi B — 3 diH p i 



d t B 
d,A 



2H, 



<n 



and fi is given by 



diHjfl + B (-3 m p i + BdiHj, 



f z = d t b--(H qi ) 2 d t B. 



(53) 
(54) 
(55) 



(56) 
(57) 

(58) 



Since we are mostly interested in the construction of flow equations, we leave the analysis 
of the physical properties of these solutions for further future investigation. 



5 Almost BPS Flow 



5.1 General Analysis 

Finally, let us investigate the almost BPS class. 

This class can be obtained from the BPS one by exploiting the very same procedure considered 
for the composite non-BPS class, the unique difference being the fact that in this case only the 
sign of the electric graviphoton charge is flipped: go - > —Qo- This yields the following first order 
system: 



diU = -Re [(3e- iao -e^ Ql ) Z t - ^s e - ll - ao+2ai) (l + e ,(OD+Jai) ) Z ti J e u , (59) 
dit = ~E l (^V3e iai Z l + y/3e i te- 2ca) Z i + 3e- iai Z ti -e iao Z i ?) e u . (60) 
By inserting (59) and (60) into the second order equations (6), a system of equations for the 
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phases «o and a.\ and for /, is obtained: 
1 



fi 



- Im 

2 



(3e~ iao +e 3iai ) Z t - V3e~ i{ao+2ai) (1 - e i{ao+3ai) ) Z tl 



-v 



<9i«o = — e Im 



e -*«i (_i + 2e 2iai + e i(ao+3Ql) ) Ili 



+ 



+ \ e U sin (oo + 3 ttl ) Wi + 1 e 2C/ cos 2 ^±^il /. 



2^3 



e 17 Im 



e -2iai (3 _ 2e 2i ai + e i(a +3ai)) jj. 



1 [7 • / , o \Txr 1 2U 2 ("0 + 3«i) 

- - e u sm (a + 3ai) Wi- - e u cos fa. 



(61) 



(62) 



(63) 



As in all other cases, fi is the imaginary part of the quantity whose real part is Wi, and it enters 
both sides of (61). 

Now, let us try to integrate out electromagnetic potential completely, as achieved in the 
composite non-BPS class. By starting from the Ansatz 



Re 



be 3U -lm 



(3. 



-iao 



e 3iai ) - > /3e- i(ao+2ai) (l + e 



»(«o+3oi)j V P 



(3e~ iao + e 3iai ) - V3e- i{ao+2ai) (l - e l(a ° +3ai) ) uf 



and inserting it into equation (61), one achieved the expression 



Im 



(3e- iQ0 + e 3iQl ) Zi - V3e- i{ - ao+2ai) (1 - e i{aa+3ai) ) Z u 



0. 



(64) 



(65) 



The equations for electromagnetic potential can be solved simultaneously iff one of these two 
conditions is fulfilled: 

I : «o + 3 oti = it, 

(66) 

II : «o + 3 oi\ = 2 arctan(6 e 2U ) . 

Correspondingly, two sub-classes I and II of full-fledged almost BPS first order equations exist. 
As we will see in the treatment below, sub-class I corresponds to BPS-like configurations with a 
positive quartic invariant of the sum of all charges pertaining to each center (namely 7 4 (^ a Q a ) > 
: BPS black hole in the single-centered limit), whereas sub-class II corresponds to non-BPS 
composite-like configurations with a negative quartic invariant of the sum of all charges pertaining 
to each center (namely I4 (%2 a Q a ) < : non-BPS black hole in the single-centered limit). 
The system of flow equations pertaining to sub-class I reads: 



diU = Re[e~ iao Zi\e u + -bf ie iU , 

dit = — e 2i ^ (e u (i + be 2U ) fi + 2 i Im [e'^Z^ ) e u E l + e iao Z u e u 
fi = d i b-2e~ u Im[e- iao Z i ] +V3be u Re\V3e- iao Zi + ie- i2 ^Z ti 



U Tpt 



(67) 
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where the phase ao satisfies the condition (already implemented in the system (67)) 



Im 



t i 



0. 



(6$ 



By setting b = in the system (67) (and using the condition (68)), one obtains the BPS first 
order system (22) constrained by (68). This leads us to identify this limit of sub-class / with 
BPS-like configurations. This conclusion is also confirmed by the asymptotical analysis, in which 
the rotation vanishes fi = 0, and b becomes zero (more rapidly than any other functions, due 
to nature of rotation); by taking (68) into account, one then gets a correspondingly constrained 
BPS system (22), consistent with I4 (^2 a Q a ) > in the single-centered interpretation; this will 
be considered in one of the examples of the next section. 

On the other hand, the system of flow equations pertaining to the sub-class // reads: 



diU 
Bit 

fi 



^Re 



e 3iai Zi - V3e iai Z-, 



t i 



-be 3U lm 
2 



-,3iai 



V3e 



n i 



X -e u E l [Vle u - 2i ^ (i - be 2U ) fi + y/3e iai Zi + e" ia 'Z t ^ + 



V3e 3iai Z { 



e iai Z 



t i 



ibe 



2U-2ia 



d;b + V3e~ 



v 



Im 



V3e 3iai Zi -e iai Z tl 



b 2 e 4C/ 
be 2U Re 



^fle^ ia Zi- e~ iai Z u 



^e 3iai Zi-e iai Z u 



(69) 



where the phase a\ satisfies the condition (already implemented in the system (69)) 



Im [e 3iai Zi] 



V3 be 2U 
2 l + b 2 e 4U 



Re 



V3 e 3iai Z, - e iai Z-, 



t i 



be 2U lm 



V3 e 3iai Z, - e iai Z-, 



t i 



(70) 



At spatial infinity, the sub-class II first order system (69)- (70) reduces to the non-BPS composite 
one (41)-(42), when fi = 0, b =constant and taking (70) into account. For this reason, as 
mentioned above we identify this sub-class as a composite non-BPS like one. 

We leave the detailed analysis of the solutions (and related constraints) of the general almost 
BPS sub-classes determined above to future investigation. The next section will be devoted to 
some illustrative examples. 



5.2 Some Examples 

1. As a first example, we consider the single-centered black hole before integrating the elec- 
tromagnetic potential out. It should be pointed out that the almost BPS class admits two 
possible single-centered limits (namely, the BPS and non-BPS one), in contrast to the BPS 
and composite non-BPS classes. By setting /j = in (59) and in (61), and replacing the 
generalized central charge and its derivative by the standard ones, one obtains 

e- = Z^f^l±ll (7 i) 
-V3Z + e 2ia ^Zt 
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for a non-BPS single-centered black hole, and 

. (e^Z + V3Z t ) 



-e iao (72) 



Z+V3e 2ia ^Z { 

for a BPS single-centered black hole. In particular, (72) is cubic in a±; its solution is 
straightforward but cumbersome, and we will not present it here. 

On the other hand, after integrating the electromagnetic fields out, the single-centered limit 
yields the following results. /) For case / of (66), the single-centered non- .BPS' limit yields 
a "small" black hole. Indeed, from / of (66), (71), (68) and from condition fi = (where 
fi is defined by (67)), the following restrictions on special geometry invariants [16] follow: 

i, = ^=(ZZ?-ZZ?)=0, (73) 

i 2 = Z t Z { = 3h = 3ZZ, (74) 
which in turn imply the quartic invariant [16] 

/4 = (*i-* 2 ) 2 -^ + 4z 4 (75) 

to vanish. On the other hand, the single-centered BPS limit yields a "large" solution; by 
combining / of (66), (72), (68), conditions fi = and 6 = 0, one can indeed check that a 
" "large" BPS black hole solution, constrained by the additional condition (73) is obtained. 

II) The same procedure applied to case // of (66) yields no "large" BPS black holes (all 
special geometry invariants are zero), whereas the non-BPS single-centered limit yields the 
same solutions as the composite non-BPS class, with no additional constraints. 

As a second and final example, we consider the well known almost BPS solution [5], [6]. 
This is obtained by constraining the phase ot\ in the almost BPS sub-class II (69) as 
follows: 

«i = — arccot [fre 2C/ ] . (76) 

The origin of such a constraint can be better understood by investigating the differential 
equations for ao and ai, and noticing that «o + «i = 7r is a particular solution. Combining 
this with the second of (66), one gets (76). By inserting this value of ati into (69) and 
considering equations for the combinations e~ u j yfy and x + bye u , the following solutions 
are obtained: 

b e~ 2U 
2 HpO Z\ 2 HpO Z\ 



e 



-4C/ _ a TJ . v3 7.2 



AH p oZf-b\ (7* 



in which the functions K\ (harmonic) and Z\ (non-harmonic) respectively satisfy the equa- 
tions 

d t K x = d * Hp ° Kl ~ dlHpl , (79) 
HpO 

d t Z x = ld l H qi +2d t K l H p oK l -d t H p oK 2 1 . (80) 
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As for fi, it reads 

f i = d i b-6H lfl Z 1 a i K 1 , (81) 
where the harmonic functions are constrained by 

diH qo = 3 diK x Zi - 3 K x + 3 K\ d i K l H p0 - diH p0 K 3 , (82) 

yielded by (70). 

We leave the analysis of the physical properties of these solutions, along with an investigation of 
new classes of solutions, for future work. 

6 Conclusion 

In the present paper, the first order formalism for multi-centered and/or rotating black holes 
was constructed and investigated in full generality for the so-called T 3 model of N = 2, D = 4 
ungauged Maxwell-Einstein supergravity. The exploited procedure sets non-BPS classes on the 
same footing as the well known BPS one. 

In order to highlight the various steps of our approach to the construction of the first order 
flow equations, we first considered the BPS class [9]. 

Similarly to the BPS class, the composite non-BPS class [3] is endowed with only one auxiliary 
phase field, which generally satisfies some algebraic constraint. For a single-centered static black 
hole, this phase can be integrated out without putting any restriction on the charges, and the 
well known non-BPS fake superpotential for the T 3 model is retrieved. It is worth stressing that 
a smooth single-centered limit is obtained also after integrating electromagnetic fields out. For 
what concerns rotating single-centered solutions, we recovered a generalization of the Rasheed- 
Larsen black hole [15], in which all electric and magnetic charges are switched on. As for 
multi-centered solutions, we also have retrieved all known solutions. 

Consistent with its nilpotent orbit characterization [3], the almost BPS class exhibits the most 
involved structure, in which two phase fields are present. We have shown that two general sub- 
classes (/ and II) exist, depending on two possible constraints satisfied by the two phases. Also, 
the almost BPS class exhibits both a BPS and a non-BPS single-centered (extremal) black hole 
limit. In particular, the BPS single-centered limit yields some restrictions on the special geometry 
invariants (explicitly derived in the discussion of some examples). This is to be contrasted with 
the BPS and composite non-BPS classes, which admit only one single-centered (extremal) black 
hole limit (BPS and non-BPS, respectively), but without any restriction on the special geometry 
invariants. Furthermore, the classes I and II of almost BPS flow respectively correspond to a 
positive and negative quartic [/-duality invariant polynomial J4, evaluated on the total charge 
vector of the system. The well known almost BPS multi-centered solution was also recovered, 
and it was shown that it is nothing but a particular solution of sub-class II. 

Having derived the most general sets of first order flow equations (consistent with the flatness 
of the three dimensional base-space), we leave the detailed study of their various new classes of 
solutions, along with their physical properties, to forthcoming works. The T 3 model provides a 
good framework, due to the presence of only one modulus and the absence of any flat direction. 
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Furthermore, we plan to extend all such results to the STU model, which instead exhibits flat 
directions [17]. 

Note added: the new paper [18] on this subject appeared on arXiv after the present work had 
been completed. 
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